Nesting symmetries and diffusion in disordered rf-wave superconductors 
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The low-energy density of states (DOS) of disordered 2D d-wave superconductors is extremely 
sensitive to details of both the disorder model and the electronic band structure. Using diagram- 
matic methods and numerical solutions of the Bogoliubov-de Gennes equations, we show that the 
physical origin of this sensitivity is the existence of a novel diffusive mode with momentum close 
to {n,Tv) which is gapless only in systems with a global nesting symmetry. We find that in generic 
situations, the DOS vanishes at the Fermi level. However, proximity to the highly symmetric case 
may nevertheless lead to observable non-monotonic behavior of the DOS in the cuprates. 



Introduction. An understanding of the quasiparticle 
(QP) excitations in the d-wave superconducting state 
of the high-Tc superconductors is essential for the elu- 
cidation of transport properties, for determining how the 
ground state deviates from the BCS model, and for de- 
scribing the instability of the lightly doped antiferromag- 
netic state to superconductivity. It has been known for 
some time that the influence of disorder on the QP states 
is quite different from ordinary superconductors, in part 
due to the gap symmetry and in part due to low dimen- 
sionality. Nersesyan et al. have shown that these two 
features conspire to introduce logarithmic divergences in 
all orders of the perturbation theory Q . Since then, sev- 
eral groups have attempted nonperturbative treatments 
of the "2D dirty d-wave problem" , arriving at a surpris- 
ingly diverse set of results. The proposed scenarios pre- 
dict vanishing ffj-^, constant HH], and divergent 
density of states (DOS) i/(e) as e ^ (energies are mea- 
sured from the Fermi level) for apparently similar models. 
Recently, two of the authors argued § on the basis of 
numerical studies that the d-wave superconductor is fun- 
damentally sensitive to "details" of disorder, as well as 
to certain symmetries of the normal state Hamiltonian. 
While this approach was successful in unifying the vari- 
ous analytical treatments, it failed to provide a physical 
explanation of the origin of this lack of robustness. 

In this work, we combine perturbative analytical and 
numerical calculations with the intent of clarifying the 
physics of the various asymptotic results for the DOS of 
noninteracting QPs. We show that the generic result is 
v{t ^ 0) — > 0, but that both the constant and divergent 
DOS can be obtained in the presence of global "nesting" 
symmetries (GNSs) which produce additional diffusive 
modes with momenta close to Q = (tt, tt) These 
symmetries drastically change the low-energy DOS and 
are explicitly obeyed by models considered in Refs. 
and where the constant and the divergent v{e — > 0) 
have been found. The GNS relations necessary for the 
TT- modes do not occur explicitly in real cuprates, but we 
demonstrate that even proximity to global nesting has 



observable consequences on thermodynamic properties. 

Models of pure d-wave superconductor. The weak- 
localization (WL) calculations will be compared directly 
with exact solutions of the Bogoliubov-de Gennes equa- 
tions on finite tight-binding lattices, details of which have 
been published elsewhere . We investigate three dif- 
ferent model band structures, each having specific sym- 
metry properties. The nearest neighbor (Nl) model with 
normal-state dispersion = — 2t(ca; -f Cy) — /i, where 
t is the hopping matrix element and is the chem- 
ical potential, is the most extensively studied. Here 
Cx.y — cos{kx,ya) and a is the lattice constant. We also 



consider a model including second-neighbor (N2) hop- 
ping, = —'2t{cx + Cy) — 4:t'cxCy — ^ aud a peculiar third- 
neighbor (N3) hopping model ^k = — 4t"(c^ -I- — 1) — /x. 
In all cases, the d-wave Hamiltonian assumes pairing be- 
ll tween nearest neighbour sites, Ak = Ao{cx — Cy). The 
lines of zeros for Ak intersect the Fermi surface ^k = 
at four symmetric nodal points for both the Nl and 
N2 models, and at eight points for the N3 model. In 
the vicinity of these points we may make an expansion 
of the QP spectrum £k — [(vFk)^ -I- (vgk)^]^/^, where 
Vi? — (9^k/c?k)k„, Vg — (9Ak/9k)k„ are the Fermi and 
"gap" velocities, and k is the momentum measured from 
the node at k„. Gapless nodal excitations determine the 
low-energy behavior of the clean DOS, Vci{e ^ 0) = ae, 
where a = N/{2TTVFUg) and N is the number of nodes. 

Self-consistent T-matrix approximation. In sharp con- 
trast with a normal metal, the DOS in d-wave supercon- 
ductors is strongly affected by disorder. The commonly 
used self-consistent T-matrix approximation (SCTMA) 
yields finite DOS at the Fermi level, — -^cijl (see, e.g., 
[pdl). Here 7 is the impurity-induced relaxation rate of 
nodal excitations, I — ln(A/7), and A ^ {vpUgY^^a^^ is 
the high-energy cut-off in the nodal expansion. 

Let us outline first the SCTMA for a general asymmet- 
ric band. Assuming that both the time-reversal and spin- 
rotational symmetries are preserved, consider a random 
distribution of point-like spinless scatterers of arbitrary 
strength U and concentration rii. The self-consistent T- 
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matrix T{e) = J2i Ti{e)Ti, depicted in Fig. |l|(a), has two 
Nambu components given by T3(e)±ro(e) = [cT5o(e)]~^j 
with c = U^^ — gz{e- = 0) and e = e — Eq. Here ti,2,3 are 
the PauH matrices and tq is the unity matrix, Sq = niT^, 
and the chemical potential shift is (5/x = S3 = riiT^. 
The quantities go, 3 are the Nambu components of the 
integrated Green's function g{e) = X]k^k(e), where 
Gk(e) = (cTo + Akfi + ekT3)/(e^ - 4)- The above, 
closed system of equations apply to both the retarded 
(R) and advanced (A) channels, and it is understood 
that e lies in the upper or lower half-plane as appro- 
priate. We remark that the quantity 33(6 = 0) is real 
and vanishes for a perfectly symmetric band; it therefore 
cannot be estimated in a nodal expansion. At the level 
of SCTMA, 33 just renormalizcs the effective scattering 
potential with U — > c^^. The Born and unitarity limits 
correspond to c ^ \go\ and c respectively. The QP 
self-energies (QPSEs) are non-singular at small frequen- 
cies e < 7 where they obtain the form S^'"^' (e) = Xe^ij. 
Here the relaxation rate 7 is the solution (cf. | |Tl| , [T^ ) of 
7 = iniTQ{e = 0), and the mass renormalization satisfies 
A = nideTo{e)\e=o- 

Analysis of diffusion modes. Since the SCTMA is ex- 
act only for a single-impurity problem, is subject to 
further corrections caused by multiple impurity scatter- 
ing which result in new low-energy regimes [0-^. We 
explore these regimes by calculating the WL correction 
to the SCTMA result in the diffusion approximation. We 
observe that the first-order correction predicts many as- 
pects of the low-energy behavior and qualitatively agrees 
with numerical calculations. In particular, the sign of the 
correction can be associated with the further tendency for 
v{e ^0): to diverge (positive sign), to vanish (negative 
sign), or to saturate (correction is absent). 

Consider the diffuson 2?q(e,e') and Cooperon Cq(e,e') 
ladder diagrams depicted in Fig. |^(b). The 2-particlc 
irreducible vertex Xi^ii — niT{e) ® T{e'). The solution 
of the matrix equations for T) and C is simplified if we 
make the decomposition Xabcd = ^J2ij ^'^''''■adTJcb 
X = T, V, C. In this basis, the ladder diagrams reduce 
to a 4 X 4 matrix equation. Furthermore, it is easy to 
show that Cq^ (e, e') = ^'q (e, e'); and that only the diago- 
nal components of the diffusive modes may be singular, 
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FIG. 1. (a) Quasiparticle self-energy in the self-consistent 
T-matrix approximation; (b) Diagramatic equation for the 
Cooperon C (the direction of the arrows in the bottom line 
should be reversed for the diffuson O); (c) diffuson and (d) 
Cooperon contributions to DOS. Solid lines are quasiparticles, 
and dashed lines indicate impurity potentials U. 

For diffusons and Cooperons, q is the difference and 
sum of momenta of QPs in the top and bottom lines in 
Fig. [|(b) respectively, and the diffusion coefficient Dq is 
found to be Dq = v^/i2-fl) with = ^{vj. + v"^). The 
logarithmic factor I in the denominator of Dq originates 
from renormalization of the real part of the QPSE. Upon 
this renormalization, Dq and vq obey the Einstein rela- 
tion Dqi^q — Ogj s^ — ■^ctv'^, where as is the universal 
value of the Drude spin conductance |||,|ll| , and s — i is 
the spin of electron. 

Due to the intrinsic particle-hole symmetry of a super- 
conductor, T2Gi(;(— e)T2 — — Gk(e), diffusive modes also 
exist in the RR and AA channels jl^ . These anomalous 
modes appear in D^^ and C^^ when e, e' < 7, 
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They are always gapless due to the fact that the spin and 
the energy are conserved quantities in the problem. 

WL correction in generic case. In the absence of spe- 
cial symmetries, Eq. (|^) gives the only diffusive modes 
which contribute to the DOS. The first-order WL correc- 
tion can be written in the form 



S^(e) = -ilm^Tr[G^(e~)E«(?)G^(g)]. 



(3) 



where ly^, = ^niTT[T{e)TiT{e')Tj]. Besides the singular 
contributions, diffusive ladders contain the regular ones. 
For superconductors, the need to sum up the nonsingular 
ladders has been pointed out in Ref. [Q. 

In the RA channel, and have diffusive poles at 
small q and e — e', 
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ttvq Dqq'^ ~ i{e - e') ' 



The skeleton diagrams for the self-energy S^(e) are de- 
picted in Figs. |^(c) and (d), with sohd grey blocks de- 
noting diffuson or Cooperon propagators. Diagram ^(c) 
- the diffuson with one leg decorated by an irreducible 
vertex, as well as a similar diagram with the other leg 
decorated - vanishes, so the leading order contribution 
is the Cooperon diagram, ^d), which yields the negative 
logarithmic correction (as found in Pp^]) 
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This holds for generic bands and arbitrary strength of 
disorder. Since DqI = iP' /(2'y), the high-energy cut-off 
A enters Eq. via the scattering rate 7 only. Further- 
more, the insertion of nonsingular ladders into diagrams 
Fig. [1(c)- (d) in all possible ways leads to ballistic renor- 
malization of the parameter Dq in Eq. (^. After the 
renormalization, Dq D, this expression is the com- 
plete singular contribution to the DOS in first order of 
perturbation theory in the inverse spin conductance crj^. 

The supression of the DOS relative to Vq at low ener- 
gies indicated by Eq. (0) is seen in nearly all published 
numerical work ]^ Jl0| , |l5| , as well in most cases presented 
in Fig. ||. A glance at Fig. ^(b)-(e) shows that, although 
j/(e 0) — !■ in these cases, there is typically more 
structure than is contained in Eq. (^). In order to un- 
derstand this better, we consider the effects of certain 
nesting symmetries. 

Global nesting symmetry T2- The above Goldstone 
modes Eq. (|^) occur for arbitrary renormalized poten- 
tial c~^. As it is seen from Eq. (Q), additional pseudo- 
Goldstone modes appear in the unitarity limit for systems 
which satisfy a GNS, 



T2GkT2 — Gk+Q. 
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When Eq. (|) is exact and c — (e.g., in the fi — 0, 
f/ = cxD Nl model studied in the propagators and 
Cq have additional poles at momenta close to Q = (tt, tt). 
These modes are related to the usual ones by 
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However, these modes are gapped by any distortion of 
the band which destroys the nesting symmetry or by the 
choice of nonunitarity potential (c 7^ 0), as we now show. 

Firstly, we note that the Nl model at /i, c = satisfies 
the GNS conditions exactly, so the 7r-modes are gapless in 
this case. Deviation from the half-filling (/i 0) creates 
a gap which can be estimated in the nodal approxima- 
tion as (5 ~ 2/1^/(7?), whereas deviation from unitarity 
(c ^ 0) yields 5 w 7(4c/7ri'o)^- Secondly, the N2 model 
has nodes at (±7r/2, ±7r/2) as n — 0, and the symme- 
try relation (j^) is approximately satisfied near the nodal 
points. However, this is not the case in the entire Bril- 
louin zone where the i'-term evidently violates Eq. (j^). 
As a result, the momentum regions located far away from 
the nodal points (these are momenta responsible for the 
global band asymmetry) contribute to a gap 5 oc t'^ in 
N2 model even at c = 0. We stress that this contribution 
cannot be estimated in the nodal approximation. There- 
fore, we conclude that the GNS relation (||) must hold 
througout the Brillouin zone, and not just near the gap 
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FIG. 2. Numerical results for m = 0.06 and Ao = 0.8. 
Solid curves are solutions of BdG equations on a 30 x 30 lat- 
tice [except for (f), which is on a 45 x 45 lattice], dashed 
curves are SCTMA results. Curves are Nl with {p, U, c} = (a) 
{0,10", 10"*}, (b) {0.6, 10*, 0.07}, (c) {0,10,0.10}; N2 with 
(d) {0, 10*, 0.099}, (e) {0, -11.675, -8 x lO"'^}, and N3 with 
(f) {0, 10,0.15}. Energies are measured in units of t for Nl 
and N2 models, and t" for the N3 model. The value of t' for 
N2 model is chosen to be 0.25. 

nodes; both deviations from Eq. (|^) and deviations from 
unitarity gap the 7r-modes. 

WL correction in ti case. The 7r-mode Cooperon 
makes a logarithmic contribution to Si>{e) of equal mag- 
nitude to Eq. (^ but - because it appears with different 
Nambu components - of opposite sign. The diffuson in 
Fig. ^c) makes a positive logarithmic contribution so the 
overall tendency is for Siy{e) to be positive. Thus, the ad- 
ditional symmetry (^) changes the asymptotic behavior 
from vanishing to divergent. If the tt- modes are gapped 
but (5 < 7, then we identify two regimes: e > (5/2 where 
gapless ordinary modes and the gapped 7r-modes con- 
tribute on equal footing such that the total correction 
increases, and e < 6/2 where the 7r-mode contribution 
saturates whereas the ordinary one continue to suppress 
the DOS. Generally, the sum of these two terms can be 
written as follows 
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Numerical calculations [p yTc| , p^ , as well as nonperturba- 
tive analytical work on the half-filled Nl model ||^ sup- 
port this result. Fig. ||(a) shows a clear peak at e = 
for the Nl model with fi,c — 0. In Figs, ^(b)-(e) the 
TT-modes are gapped in three different ways. In all cases, 
the qualitative behavior is the same: there is an initial 
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upturn starting at e ~ 7, which then crosses over to a 
downturn. As we move farther from the symmetric uni- 
tarity Hmit, either by increasing ji, decreasing U or in- 
creasing second-neighbour hopping, the peaks are shifted 
to higher energy, then become less pronounced, and af- 
ter that disappear. In (b)-(e), the shape of the curves 
is qualitatively what one might expect based on Eq. (|^). 
While this expression regularly overestimates the posi- 
tion of the maxima predicted as e ~ (5/4 (cf. Fig. |^), 
it is clear that the qualitative explanation based on the 
first-order perturbation correction could not pretend to 
describe all the details of the numerical experiment. We 
mention finally that the numerics [see Fig. |^(e)] and our 
WL calculation are in agreement that it is not possible to 
produce a divergence by means of fine-tuning parameters 
in the non-nested case, as suggested in Ref. [t). 

Models with other nesting symmetries. There is a sec- 
ond possible nesting condition. 



TsGkTa = Gk^ 
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which leads to a non-generic DOS. A model (N3) with 
this symmetry which yields an exact solution was recently 
discussed in Q . In this instance it was found, at least for 
Lorentzian disorder, that i/(e) remains finite as e — > 0. A 
similar result is seen in the WL calculations. While the 
0- modes are unaffected by Eq. (0), the 7r-modes change 
drastically. Divergences now appear in D^^ and with 
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In this case, the diffusive modes contribute for arbitrary 
c. Calculating the WL correction to the DOS we ob- 
serve, as before, the exact cancellation between 0- and 
TT-Cooperons. However, the contribution of vr-diffusons 
[Fig. |^(c)] now also sums to zero, and the total first-order 
WL correction to DOS is Si'{e) = 0. 

The N3 model explicitly satisfies the symmetry (|^) for 
arbitrary /i, and a typical case is shown in Fig. ||(f). As 
expected, the DOS is finite at e = 0, and agrees closely 
with the SCTMA result. This resembles the situation in 
a noninteracting normal metal [ 3|, where the disorder by 
itself does not affect the DOS [1^. On the other hand, 
this model exhibits a Meissner effect. It is clear, how- 
ever, that Eq. describes a specific situation which is 
probably unrelated to the real cuprates [^ . We note that 
the approach of Ref. [Q fails to predict the characteristic 
behaviors of the f2 and models identified here. 

In order to complete our analysis of various extra sym- 
metry possibilities, let us mention the purely artificial 
case: riGk^i = Gk+q- Then the gaplcss tt- modes appear 
in the unitarity regime only, for which case we predict the 
constant ^{e 0), since 0- and 7r-Cooperons cancel each 
other whereas the diffuson contributions sum to zero. 

In conclusion, we investigated the low-energy den- 
sity of states in a disordered 2D d-wave superconductor 



and argued that the generic feature of DOS is its low- 
energy suppression. However, due to proximity to perfect 
nesting, pseudo-Goldstone diffusive mode with momenta 
close to (tt, tt) gives rise to an energy dependence of DOS 
which can be strongly non-monotonic. The physics of 
this novel mode appears to account for the differences 
between many recent nonperturbative treatments and ex- 
plains numerical results. The effects we predict should be 
visible in low-temperature specific heat experiments on 
disordered cuprates. In layered compounds with weak 
interlayer coupling t±, WL corrections to the DOS ap- 
pear for t± < 7, but saturate for e < t'^/7. Thus the 
criterion for a pronounced energy window in which WL 
corrections occur is i^/7 < e < 7. This regime should 
be accessible in underdoped highly disordered materials. 
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